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Abstract: Since the discovery of the halo nucleus !!Li in 1985, halo phenomena in exotic nuclei have always been 
an important frontier in nuclear physics research. The relativistic density functional theory has achieved great success 
in the study of halo nuclei, e.g., the self-consistent description of halo nucleus !!Li and the microscopic prediction of 
deformed halo nuclei. This paper introduces some recent progresses, including the investigation of halo nucleus "Mg 
and the prediction of the N = 28 shell collapse and a deformed halo in the new isotope ??Na based on the deformed 


relativistic Hartree-Bogoliubov theory in continuum (DRHBc), as well as the exploration of triaxially deformed halo 


nuclei by the newly developed triaxial relativistic Hartree-Bogoliubov theory in continuum (TRHBc). 
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Introduction 


The discovery of the first halo nucleus"! in 1985 has 
attracted great interest. The radius of a halo nucleus sig- 
nificantly deviates from the empirical formula rA" 3 (where 
rg % 1.2 fm and A is the mass number), shaking the sim- 
ple assumption that the nucleus is incompressible. Therefore, 
the description of halo phenomena has challenged traditional 
models for nuclear structure. The experimental exploration 
of halo nuclei also has been stimulating the development of 
new generations of radioactive ion beam facilities worldwide. 

During the past decades, about 20 halo nuclei or can- 
didates have been identified or suggested in experiments, as 
shown in Fig. 1. On the theoretical side, the nuclear models 
studying halo phenomena include the few-body model [%3], 


shell model'**!, antisymmetrized molecular dynamics ©”), 


[8-9] [10-11] 


, nonrelativistic and rel- 


[12-13] 


halo effective field theory 
ativistic density functional theories , etc. 


The relativistic density functional theory is designed for 
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nuclear many-body problems based on the quantum field the- 


[15] 


ory and density functional theory ™. Due to many appeal- 


ing advantages €, such as the automatic inclusion of the 


spin-orbital coupling 5-1?! 


and the reasonable description of 
its isospin dependence"?l, the competition between scalar 
and vector densities giving rise to the new saturation mecha- 
nism", natural explanation of the pseudospin symmetry in 
the nucleon spectrum ???^! and the spin symmetry in antin- 


[24-26] and self-consistent treatment of the 


[27-28] 


ucleon spectrum 
nuclear magnetism , the relativistic density functional 


theory has become one of the most popular nuclear theo- 
ries [1529-38]. 


In the study of halo nuclei, great progress has been 
achieved based on the relativistic density functional theory. 
The relativistic continuum Hartree-Bogoliubov (RCHB) the- 


ory [12,39] 


, which solves the relativistic Hartree-Bogoliubov 
(RHB) equation in coordinate space by assuming the spheri- 
cal symmetry, is suitable for describing spherical exotic nu- 
clei. Based on the RCHB theory, the halo phenomenon 
in !!Li was described in a microscopic and self-consistent 
way!!?], the giant halo formed by more than two neutrons 
was predicted", and the neutron halo in hypernuclei was 
explored #1, 

The deformed relativistic Hartree-Bogoliubov theory in 


continuum (DRHBc) [42-45] | Which considers the axial defor- 


bir 
cn 
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Fig. 1 


(Color online) Experimentally known nuclear landscape from helium to phosphorus, where stable nuclei and experimentally con- 


firmed/suggested neutron as well as proton halo nuclei/candidates are indicated by grey, olive/green, and orange/yellow colors, 


respectively. Taken from Ref."*! and slightly modified. 


mation and solves the RHB equation in a Dirac Woods-Saxon 
(DWS) basis'“°!, can provide an adequate description for ax- 
ially deformed exotic nuclei. The DRHBc theory has pre- 
dicted the deformed halo nuclei ?^^Mg and the shape de- 


coupling between the core and the halo ?-?! 


, which is con- 
sidered one of the interesting new phenomena near the drip 
line, While the heaviest magnesium isotope synthesized 
in the laboratory so far is Mg?! and the DRHBc predic- 
tion for ^^ Mg remains to be verified, many known halo 
nuclei have been successfully described using the DRHBc 


17195 [49-50] 15.19.2512]. 31 (53) ang 


theory, including 
37Mg*951_ A deformed two-neutron halo in ?Na has also 
been predicted by the DRHBc theory!^!. In addition, a 
DRHBc mass table including both deformation and contin- 


uum effects is under construction 5659]. 


During this con- 
struction, many interesting phenomena have been investi- 
160-69], The 


DRHBc theory has also been extended with the angular mo- 


gated or predicted using the DRHBc theory 


mentum projection, enabling the exploration of rotational ex- 


citations of deformed halo nuclei "07", 


It should also be mentioned that there are many stud- 


ies on deformed halos using other models, e.g., the single- 


particle model ?!, antisymmetrized molecular dynamics ?!, 


[74-75] 


Hartree-Fock theory , Hartree-Fock-Bogoliubov the- 


[76-77] [78-79] 


ory , relativistic Hartree-Fock-Bogoliubov theory 
and ab initio theory ®°!, Most of these studies have assumed 
axial symmetry for nuclear deformation. 

The triaxial relativistic Hartree-Bogoliubov theory in 
continuum (TRHBc) [81] which considers the triaxial defor- 
mation and solves the RHB equation in the DWS basis, has 
been developed very recently in order to explore triaxially 
The observed heaviest odd-odd alu- 


minum isotope, *7A1'*!, has been predicted as a triaxial halo 


deformed halo nuclei. 


nucleus with a novel shape decoupling between the core and 
the halo at the triaxial level ®!!, 

This paper is structured as follows. The frameworks of 
the DRHBc and TRHBc theories are briefly presented in Sec- 
tion 1. Several of the above mentioned studies on halo nuclei 
are introduced in Section 2, including 37Mg， 3°Na, and ?Al. 


A summary is given in Section 3. 


1 The DRHBc and TRHBc theories 


In the following we briefly introduce the main formal- 
ism for the DRHBc and TRHBc theories. Starting from a 
meson-exchange or point-coupling Lagrangian density, the 


RHB equation, 


hp-A å UN p [Vr (i 
-4 nn +a} AV, NYS? 


J 原子 核 物 
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can be derived by treating the nuclear mean field and pairing 
correlations on the same footing *?!. In Eq. (1), A is the Fermi 
energy, and E, and (U,, V,)' are respectively the quasiparti- 
cle energy and the quasiparticle wave function. hp denotes 


the Dirac Hamiltonian, 
hp(r) 2 a-p-Ví(r)- B[M + S(r)], (2) 


with the scalar potential S(r) and the vector potential V(r). 


A represents the pairing potential, 
A(r,,r,;) = VPP(ri, r,)k(r,.r;), (3) 


where VPP is the pairing interaction and x is the pairing 
tensor. The present DRHBc and TRHBc theories adopt a 


density-dependent zero-range pairing interaction, 


pir) 


V(r, ra) = E — P^)ó(r, — r2) [ = | , (4) 


sat 
in which Vp is the pairing strength, za — P^)is the projection 
operator for the spin-zero component, and p,,, is the satura- 
tion density of nuclear matter. 

In the DRHBc theory for axially deformed nuclei, the 


densities and potentials are expanded as follows, 
f(r) = Y fi) Picos 0), (5) 
1 


where P, is the /-order Legendre polynomial, and the as- 
sumed spatial reflection symmetry limits / to being even 
numbers. 

In the TRHBc theory for triaxially deformed nuclei, the 
densities and potentials are expanded in terms of spherical 


harmonics, 


SO = 2 fin( Yin(0, o). (6) 
im 


The spatial reflection symmetry and the mirror symmetries 
with respect to xy, xz, and yz planes lead to the limitations 
that / and m are even numbers and fj,(7) = f, ,(r) 51. 

The RHB equations are solved in the DWS basis 6941, 
whose wave function has a more appropriate asymptotic be- 
havior compared with the harmonic oscillator basis and is 
thus suitable for expanding the wave functions of weakly 
bound nuclei. 

In practical calculations, an energy cutoff of Eu = 300 
MeV and an angular momentum cutoff of J max = 23/2 h for 
the DWS basis can guarantee the convergence of numerical 
results 6l. 
—325 MeV-fm?, the saturation density p = 0.152 fm^?, 


and a pairing window of 100 MeV can reproduce the odd- 


In such a basis space, the pairing strength V, = 


even mass differences of calcium and lead isotopes®®!, For 


light nuclei, it is enough to truncate the Legendre or spherical 
harmonic expansion at lax = 6/?!, The blocking effects in 
odd-mass and odd-odd nuclei are taken into account via the 


equal filling approximation ^^55561. 


2 Recent progress 


2.1 The halo nucleus Y Mg 


The deformed p-wave neutron halo nucleus 37Mg, con- 
firmed experimentally in 2014/5755]. is the heaviest nuclear 
halo system to date. Although ?"Mg had been synthesized 
and identified in 1996 [89]. theoretical studies over the subse- 
quent two decades, including phenomenological calculations 


based on a Woods-Saxon potentia] P?! 


and microscopic ones 
based on nonrelativistic and relativistic density functional 
theories ?'-?5!, did not predict its neutron halo. Even after 
2014, a fully microscopic description of the neutron halo in 
37Mg remained a challenge. For instance, a large quadrupole 
deformation (f, 2 0.5) must be introduced in the calculations 
based on Woods-Saxon potentials so as to reproduce the p- 
wave occupation of the valence neutron in "Mg^?7l. The 
calculations based on the antisymmetrized molecular dynam- 
ics reproduced the matter radii of ?^?9Mg but significantly 
underestimated the matter radius of "Mg ???!, The results 
from nonrelativistic Hartree-Fock-Bogoliubov calculations 
were found to be density-functional dependent 67799]. with 
only the M3Y-P6 functional capable of providing the p-wave 


components of the valence neutron in Mg 8l. 


The deformed p-wave neutron halo in "Mg has been 
described in a microscopic and self-consistent way using the 
DRHBc theory 4). Figure 2 illustrates the two-dimensional 
neutron density distributions for 537Mg in the DRHBc cal- 


H00] "The neu- 


culations with the density functional PC-F1 
tron density distribution of "Mg is remarkably more diffuse, 
which is consistent with the picture of nuclear halos. The 
root-mean-square (rms) radii for ?-?" Mg calculated from nu- 
clear densities are respectively 3.44, 3.49, and 3.57 fm, in 
satisfactory agreement with the empirical radii of 3.44(03), 
3.49(01), and 3.62(03) fm deduced from the reaction cross 


1731, Furthermore, the densities and the wave func- 


sections 
tion of the valence neutron from the DRHBc theory can be 
fed as microscopic inputs into the Glauber model to study 
reaction observables, which turns out successful in describ- 
ing the halo nucleus *!Ne?!, Such a work for magnesium 


isotopes has also been performed recently, and the reaction 
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Fig. 2 (Color online) Two-dimensional neutron density distributions for 3 -37Mg as well as "Mg's core and halo in the xz plane with z 


being the symmetry axis. Taken from Ref. >“). 


cross section of "Mg bombarding a carbon target and the 
longitudinal momentum distribution of Mg fragments after 


its breakup are reproduced very well 55, 


To further understand the halo phenomenon in 37Mg， 
Fig. 3 presents the single-neutron orbitals around the Fermi 
energy and their rms radii. It is found that the rms radius of 
the weakly bound orbital 5 is approximately 5.6 fm, signifi- 
cantly larger than others. Orbital 5 is occupied by the valence 
neutron in "Mg, with 38% 2p1,, and 19% 2p,,, components, 
which are comparable to the ~ 40% p-wave halo components 
suggested in experiments 7!, The low centrifugal barrier for 
p waves allows considerable tunneling of the wave function 
into the classically forbidden region, which, together with 
the weak binding of orbital 5, results in its largest rms ra- 
dius. While certain s- and p-wave components are found for 
orbitals 1 and 2, respectively, the diffuseness of their wave 
functions is suppressed by the relatively deep binding. For 
orbitals 3 and 4 dominated respectively by d and f waves and 
orbitals 6 and 7 with the occupation probabilities < 0.1, their 
contribution to the neutron halo in "Mg is marginal. There- 
fore, it is natural to decouple the neutron density of "Mg 
into a core and a halo, with the Fermi energy as a threshold. 
As seen in Fig. 2, the core density closely resembles that of 
36Mg, while the halo density predominantly contributes to 


the total neutron density far away from the center. 


The results from DRHBc calculations based on other 
density functionals, such as PC-PK1 0! NL3* 021. and 
PK1U9?l are quantitatively consistent with the above re- 
sults. Therefore, the microscopic self-consistent description 
of the DRHBc theory on the halo nucleus *’Mg is essen- 
tially density-functional independent. It is regrettable that 
the DRHBc theory was only applied to the even-even nuclei 
in magnesium isotopes before 2014, missing the prediction 


of the neutron halo in "Mg. 
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Fig.3 (Color online) Energy versus rms radius for single-neutron 


orbitals around the Fermi energy 4, (dashed line) in the 
canonical basis for "Mg. The order i and quantum num- 


bers 2" are given on the right. Taken from Ref. ^l. 
2.2 The new isotope 99Na 


The exploration of the boundaries in nuclear chart, i.e., 
drip-line locations, has always been a major goal of nuclear 
physics!??l, Although the proton drip line has been experi- 
mentally delineated up to neptunium with Z = 93"), the 
neutron drip line is only known up to neon with Z = 10061. 
In 2022, further experimental exploration discovered the cur- 
rently most neutron-rich isotope of sodium with Z = 11, 
??Na, but the neutron drip line has not been determined 
yet! ?7l. The neutron number of ?Na is N = 28, a tradi- 
tional magic number. The disappearance of traditional magic 
numbers and the emergence of new ones are novel phenom- 
ena in exotic nuclei near drip lines, e.g., the disappearance 
of the N = 28 magicity in the isotones of ?Na, '^Mg and 


4， 原子 核 物 理 评 论 


第 卷 


-4 Tj, 
o m I 
= 
E js 3/2- 
D OF 2pio T MM 
7 uec c I se 
5 Ihre 1/2 v2 
z -4t 3/2 1.0 
c N -20 12 0.8 
2 0.6 
D 
= aa 
DD -8 上 d 

1dy; 1/2* 02 

0.0 
-0.2 0.0 0.2 0.4 0.6 


Quadrupole deformation B; 


Fig.4 (Color online) Single-neutron levels around the Fermi en- 
ergy A, (dashed line) of ?Na in the canonical basis from 
constrained calculations. Their quantum numbers n/j in 
the spherical limit and 42^ on the prolate side are labeled. 
The occupation probability v? is scaled by colors. The grey 
vertical line corresponds to the ground state (g.s.) of ?Na. 
Taken from Ref. !4!. 


r,/ fm 


Fig. 5 (Color online) Neutron density distributions (a) along (0 — 
0°) and (b) perpendicular to (9 = 90°) the symmetry axis 
for neutron-rich odd-even sodium isotopes 3.3…“ Na, In 


(b), r, = Vx? + y?. Taken from Ref. U^. 


42.5; 108-109] On the other hand, as shown in Fig. 1, clues to 


halo nuclei or candidates are found in every isotopic chain 
from He to P, except for Na and Si. Therefore, it is necessary 
to study the shell structure of ?Na and explore possible halo 
phenomena in neutron-rich Na isotopes, which could provide 


references for future experiments. 


Figure 4 shows the evolution of single-neutron levels 
around the Fermi energy with the quadrupole deformation 
obtained from DRHBc calculations with PC-PK1. In the 
spherical limit, the orbitals 2p," and 1 f7,, are nearly degen- 
erate and close to the 2p,5. In the traditional shell model, 
there is a sizable energy gap between 1 f}; and 2p3,, form- 
ing the N — 28 shell closure and making the spherical shape 
energetically favored. In ?Na, the lowering of 2p orbitals re- 
sults in the collapse of the N — 28 shell closure, and the 
down-sloping levels split from the 1 f,,. orbital leads to a 
large ground-state deformation of f, > 0.4. In addition, as 
shown in Fig. 4, the pairing correlations lead to the partial 
occupation of the weakly bound or continuum 1/27 and 3/27 
orbitals by the valence neutrons in ?Na. Due to the deforma- 
tion effect, the pf components are mixed in the wave func- 
tions of these orbitals, and the p-wave components in valence 
neutrons favor the formation of a neutron halo. Therefore, the 


halo phenomenon may exist in ??Na. 


Figure 5 shows the neutron density distributions along 
and perpendicular to the symmetry axis for the odd-mass Na 
isotopes with N > 20. The ground-state deformation of *! Na 
is spherical due to the N — 20 shell closure, and other more 
neutron-rich Na isotopes are well deformed with quadrupole 
deformation fj, > 0.35. Therefore, the significant increase 
in density distribution along the symmetry axis from ?'Na 
to Na in Fig. 5(a) can be understood from the deforma- 
tion effects. From ??Na to ^'Na, the neutron density distri- 
bution along the symmetry axis gradually becomes more dif- 
fuse with the increasing neutron number. In Fig. 5(b), sim- 
ilar gradual growth is seen in the neutron density distribu- 
tion perpendicular to the symmetry axis from 31Na to 37Na. 
However, far away from the center, the density distributions 
of ?Na and *! Na are very diffuse, even though they are pro- 
lately deformed. This suggests oblate neutron halos in ?Na 
and *!Na, with the prolate-oblate shape decoupling for the 


core and halo. 


Further analysis reveals that the 1/27 and 3/2. orbitals 
embedded in continuum contribute predominantly to the ha- 


los in ?^'Na, These two orbitals have considerable p-wave 
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components. According to the angular momentum cou- 
pling, both spherical harmonic functions |Yi0(0, p)|? and 
|Y144 (0, q)l contribute to the 1/27 orbital, and the latter 
dominates. For the 3/2 orbital, only |Y;,,(6, q)l con- 
tributes. The angular distribution is prolate for |Y,)(9, Q9) « 
cos? 0 and is oblate for N41, Oo)|? « sin? 0. The shape of 
the halo is thus oblate. 

Therefore, based on the DRHBc calculations with den- 
sity functional PC-PK1, the disappearance of the traditional 
magic number N = 28 and the appearance of a deformed 
neutron halo in the new isotope ??Na are predicted. For the 
results from other density functionals for Na isotopes, the dis- 
cussion on the neutron drip-line location, and further explo- 
ration on the microscopic mechanism for the formation of 


halos, see Ref. 4]. 


2.3 The possible triaxial halo nucleus 42 Al 


Non-axial deformation, namely, triaxial deformation, is 
one of the basic deformation degrees of freedom in atomic 
nuclei. The triaxiality is closely related to rich nuclear phe- 
nomena, e.g., the chiral rotation! 10, wobbling motion HH, 
and fission !!!?!, It is recently found that the information on 
nuclear deformation can be extracted from relativistic heavy- 


i 13]. and the evidence for the triax- 


ion collision experiments 
ial structure in ?Xe has been reported by analyzing the data 
from the CERN Large Hadron Collider! ^l, which further 


stimulates the interest in triaxial nuclei. 


It would be interesting to explore the possible existence 
of triaxial halo nuclei. In 2021, a calculation based on the 
Woods-Saxon potential suggests that the region of halo nu- 
clei could be extended by triaxial deformation that allows 
the appearance of s- or p-wave components in some weakly 


[15]. However, the triaxial deformation, the 


bound orbitals 
depth, and the width of the potential are adjustable param- 
eters, and the crucial pairing and continuum effects are not 
included in such a phenomenological model. In contrast, the 
TRHBc theory takes into account the triaxial deformation, 
pairing correlations, and continuum effects in a microscopic 
and self-consistent way, making it capable of describing tri- 
axial halo nuclei properly. The TRHBc theory predicts ^AI 
as a triaxial halo nucleus by examining the neutron separation 
energies, rms radii, density distributions, and single-neutron 


orbitals around the Fermi energy for aluminum isotopes ?!!, 


Figure 6 exhibits the density distributions in xy, xz, and 


yz planes for the core and the halo of 42A], calculated from 


the TRHBc theory with PC-PK1. Itcan be found that the halo 
density is significantly more diffuse than the core density, es- 
pecially in the yz plane. Quantitatively, the rms radius is 3.85 
fm for the core and 5.26 fm for the halo. The deformation 
parameters (f, y) are (0.38, 50°) and (0.79, —23°) for the core 
and the halo, respectively. The negative y value means that y 
is the intermediate axis while x is the short one. This is just 
the reverse of the case for the core. With the corresponding 
rms radius and deformation parameters, the schematic pic- 
tures for the core and the halo are also displayed in Fig. 6, 
where the long, intermediate, and short axes can be distin- 
guished clearly. Therefore, ?AI is a triaxial halo nucleus 
with a novel triaxial shape decoupling between the core and 
the halo. For detailed results and discussion on the triaxial 


halo characteristics of 2 AT, see Ref, ?!1. 


3 Summary 


The halo phenomena in exotic nuclei have long been 
an important frontier in nuclear physics research. The rel- 
ativistic density functional theory has made great progress 
in the study of halo nuclei. Recently, a microscopic and 
self-consistent description of "Mg, including its small one- 
neutron separation energy, large rms radius, diffuse density 
distribution, and p-wave components for the halo neutron, 
has been obtained using the DRHBc theory P^... Based on 
the structure input purely from the DRHBc theory, the re- 
action cross section of ?"Mg on a carbon target and the mo- 
mentum distribution of Mg fragments after its breakup are 
reproduced very well by the Glauber model P?!, Therefore, a 
unified description of the halo characteristics of "Mg from 
nuclear structure to reaction dynamics has been achieved. 

The recently discovered isotope ?Na has been inves- 
tigated with the DRHBc theory'?!. It is revealed that the 
lowering of 2p orbitals in the spherical limit results in the 
disappearance of the N — 28 magicity in ?Na. The pairing 
correlations and the pf component mixing driven by defor- 
mation lead to the partial occupation of the weakly bound and 
continuum orbitals with certain p-wave components, giving 
rise to the formation of an oblate neutron halo around the 
prolate core in ?Na. 

To explore triaxial halo nuclei, the TRHBc theory has 
been developed, which self-consistently includes the triaxial- 
ity, pairing correlations, and continuum effects *!. The ob- 
served heaviest odd-odd aluminum isotope, ^^ AI, is predicted 


as a triaxially deformed halo nucleus, in which a novel triax- 
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Fig. 6 (Color online) Neutron density distributions in xy, xz, and yz planes for the core (a)- (c) and the halo (d)-(f) in the predicted triaxial 


halo nucleus ? Al. In each plot, a circle in dotted line is drawn to guide the eye. With the rms radius and deformation parameters 


f and y calculated from densities, the corresponding schematic shapes for the core and the halo are given in the left, in which s, i, 


and | respectively represent the short, intermediate, and long axes. Taken from Ref. *! and slightly modified. 


ial shape decoupling, i.e., the exchange of the intermediate 
and short axes between the triaxial core and the triaxial halo, 


is found. 
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摘要 : 自 1985 年 实验 发 现 晕 核 U Li 以 来 ， 奇 特 原子 核 中 的 坚 现 象 一 直 是 核 物 理 研 究 的 重要 前 沿 。 相 对 论 密度 
泛 函 理论 在 尝 核 研究 中 取得 了 很 大 成 功 ， 如 坚 核 0Li 的 自 洽 描 述 和 形变 晕 核 的 微观 预言 。 本 文 介绍 一 些 近 期 
和 进展， 包括: 形变 相对 论 连 续 谱 Hartree-Bogoliubov HH HER Mg、 预言 新 核 素 ?2Na 的 N = 28 384 
和 形变 晕 ， 以 及 新 发 展 的 三 轴 相 对 论 连 续 谱 Hartree-Bogoliubov 3836 36 R — A W 3c A. 
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